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Highly compressed nuclear matter created in relativistic 
heavy collisions is to large extent governed by local non- 
equilibrium. As an idealised scenario colliding nuclear matter 
configurations are studied with an effective in-medium inter- 
action based on the microscopic DBHF model. It is found that 
on top of the repulsive momentum dependence of the nuclear 
forces kinetic non-equilibrium leads to an effective softening 
of the equation of state as compared to ground state nuclear 
matter. The separation of phase space which is the basic 
feature of such anisotropic configurations has thereby a sim- 
ilar influence as the introduction of a virtual new degree of 
freedom. 



One major goal of relativistic heavy ion physics is to 
explore the behaviour of the nuclear equation-of-state 
(EOS) far away from saturation, i.e. at high densities 
and non-zero temperature. Over the last three decades a 
large variety of observables has been investigated both, 
from the experimental and theoretical side motivated by 
the search for the nuclear EOS. The collective particle 
flow is thereby intimately connected to the dynamics dur- 
ing the compressed high density phase of such reactions 
10,^. E.g., the elliptic flow which develops in the early 
compression phase is thought to be a suitable observ- 
able to extract information on the EOS But also the 
production of strange particles is a good probe to study 
dense matter Q]. Recent precision measurements of the 
_ftr+ production at SIS energies strongly support the sce- 
nario of a relatively soft EOS |,§. 

The temporal evolution of the collision from a highly 
anisotropic initial configuration in phase space to an - 
at least partially - equilibrated final configuration is suc- 
cessfully described by microscopic transport models like 
BUU or QMD |^ . In this type of models the nuclear 
mean field is usually based an phenomenological parame- 
terisations [^-10|. Such parameterisations allow different 
extrapolations to high densities, subsummized by refer- 
ring to a 'hard' or a 'soft' equation-of-state jl], which 
can be tested in heavy ion collisions. A fundamental 
question is, however, if such a procedure is sufficient to 
to extract well defined information on the EOS of equi- 
librated nuclear matter (NM) from heavy ion collisions. 
The mean field used in transport simulations refers to 
locally equilibrated nuclear matter. Anisotropy effects of 
the local phase space are usually neglected for the den- 
sity dependent (or local) part of the mean field although 
it is clear that they should be included from a theoret- 



ical point of view |ll|-|14|. Such a treatment appears to 
be justified if non-equilibrium effects are small, have a 
short lifetime, or if the mean field is not much affected. 
As we will discuss in the present work none of this three 
conditions is fulfilled in relativistic heavy ion reactions 
in the SIS energy range: 1. The initial phase space dis- 
tribution in the participant zone is that of two currents 
of nuclear matter colliding with beam velocity. The lo- 
cal momentum distributions of colliding nuclear matter 
configurations, called CNM in the following, are given by 
two Fermi ellipsoids, i.e. two boosted Fermi spheres sep- 
arated by a relative velocity ||l^. 2. The relaxation time 
needed by the system to equilibrate coincides more or 
less with the high density phase of the reaction. Hence, 
anisotropy effects are present all over the compression 
phase where one essentially intends to study the EOS 
at high densities. Experimental evidence for incomplete 
equilibration even in central collisions at SIS energies has 
recently been reported in |l^. 3. The impact of phase 
space anisotropics on the nuclear EOS is large and leads 
to a considerable softening of the effective EOS seen in 
heavy ion collisions. 

To obtain a more quantitative measure for the size and 
relevant time scales for phase space anisotropics in Fig. 
1 the time evolution of the quadrupole moment of the 
energy-momentum-tensor 
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at the collision centre of central (b=0 fm) Au + Au re- 
actions at two typical beam energies, i.e. 0.6 A.GeV and 
6.0 A.GeV, respectively, is compared to the time evolu- 
tion of the corresponding baryon density. The quantities 
were obtained in relativistic BUU calculations jlj]. Qzz 
is a measure for the anisotropy in beam (z) direction. At 
6.0 A.GeV Qzz reaches already the asymptotic value of 2 
where the Fermi momentum can be neglected compared 
to the beam velocity. Fig. 1 demonstrates that the relax- 
ation time to reach equilibrium configurations {Qzz — 0) 
coincides more or less with the high density phase, inde- 
pendent of the beam energy. Therefore non-equilibrium 
effects should be taken into account on the level of the ef- 
fective in-medium interaction which means to determine 
the mean field used in transport calculations consistently 
for colliding nuclear matter. 
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FIG. 1. Time evolution of the quadrupole moment Q^^ 
(bold lines) at the collision centre in central Au-\-Au reactions 
at 0.6 and 6.0 A.GeV. Qzz is a measure for the anisotropy of 
the local momentum space during the reaction. For compar- 
ison also the central densities (in units of psat) are shown. 

In the present work we want to study the impact of 
such non-equilibrium effects on the nuclear EOS probed 
by heavy ion collisions. The configurations are idealised 
by colliding nuclear matter which describes well the lo- 
cal phase space in the participant zone of the collisions. 
CNM also provides a smooth transition to equilibrated 
configurations which evolve in the later stages of the re- 
actions . The effective interaction used in the present 
studies is based on the microscopic relativistic Dirac- 
Brueckner-Hartree-Fock (DBHF) approach [ pT[p^ which 
turned out to be quite successful in the description of nu- 
clear matter bulk properties. Here the nucleon-nucleon 
interaction is fixed by free NN-scattering and the medium 
correlations contained in the self-consistent summation of 
the ladder diagrams of the Bcthe-Salpeter equation de- 
termine the behaviour of the nuclear EOS. DBHF forces 
have been extensively tested at SIS energies below 1 
A.GeV, and a general agreement with corresponding flow 
data has been observed The mean field used 

in the present work is based on recent DBHF calculations 
psf with nuclear matter saturation properties (Bonn A 
potential) of p^at = 0.185 fm-^ S'''"^ = -16.15 MeV 
and a compression modulus oi K = 230 MeV. The so- 
lution of the full DBHF problem for colliding nuclear 
matter, i.e. two-Fermi-ellipsoid configurations, is, how- 
ever, still an unresolved problem and thus we determine 
the effective interaction in CNM as discussed in [ p^ . 
By this procedure the density and momentum depen- 
dence of the DBHF self-energies are extrapolated to the 
CNM configurations and the mean field is constructed 
by a superposition of the contributions from the two cur- 
rents. However, already this approximation invokes a 
self-consistency problem for the CNM configuration 

912 = 61 + 02-91 -62 (2) 

composed by the two currents 9^ = Q{ii*{kp.) — k*u^). 
The effective mass M* = M + £5^3 , the scalar density 
gsj^2 =< M*/E* >i2 and the configuration (||) are cou- 
pled by non-linear equations. 9 is the step function. 



kp. are the Fermi momenta and itj' = (71,71 u^) are the 
streaming velocities of the two subsystem currents. The 
last term in eq. (||) ensures that the Pauli principle is ful- 
filled for small velocities where the two ellipsoids might 
overlap. For details see Jl3t . 

< X >i2= J d^kX(fc)9i2(fc) (3) 

denotes the summation over all occupied states. In spin- 
isospin saturated nuclear matter the phase space occu- 
pancy factor is K = 4. The energy momentum tensor in 
CNM is given by [|l| 

Tf^"^ =< fc'^rVi?* >i2 ^V.r - Qs,, - V>^,,} 

(4) 

with the scalar contribution 

-< ^s,,M*/E* >i2 /qs,, (5) 

and the terms arising from the vector field 

=< ^'i2k*''/E* >i2 . (6) 

The scalar self-energy £5^3 enters in a Hartree form 
averaged over the explicit momentum dependence into 
(^. In an analogous way a Hartree vector self-energy 
^12 — ^12^ ii2\/ ji2 can be defined. The full self-energy 
has the form Si2(fci?j , /ci?^ , ui, M2) = '^Svi + 7ai^i2- 
contrast to the local density approximation where the 
self-energy is a function of the total Fermi momentum 
kpt^t the self-energy in CNM depends on the densities of 
the two subsystems and their streaming velocities [|l3|,0 . 

In the following we will only consider the symmetric 
case (fc^i — kp^)- In Fig. 2 the equation-of-state's in 
symmetric colliding nuclear matter are shown for stream- 
ing velocities of the subsystem currents u = |ui| = 
0.2/0.4/0.6/0.8 (in units of c) which correspond to in- 
cident laboratory energies of i?iab=0. 08/0. 36/1. 05/3. 34 
GeV/nucleon. u = is the isotropic case (ground state 
NM). The energy per particle is defined in the usual way 
as 

Si2(to,") = T™/ei2-M (7) 

with Q12 = \/ ji2 1 >12 \c.m. the invariant baryon 
density in the cm. frame of the two currents, i.e. the 
frame where ji2 — 0. At high streaming velocities the 
"EOS" is significantly stiffer then in ground state nuclear 
matter because the energy per particle £^12 includes the 
contribution form the relative motion of the two currents. 
The general repulsive character of the momentum depen- 
dent part of the nuclear interaction leads to a strongly 
enhanced repulsive component which increases with an 
increasing amount of relative motion in the system. 
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FIG. 2. EOS in nuclear matter (solid) and colliding nu- 
clear matter determined in the DBHF model. The upper 
part shows the total energy per particle E12 as a function of 
the cm. total density. The streaming velocities are u—0.2 
(dotted), 0.4 (dashed), 0.6 (long-dashed), 0.8 (dot-dashed). 
The lower part shows the effective EOS, i.e. the binding en- 
ergy per particle -Ei"^"* where the kinetic energy of the relative 
motion in CNM has been subtracted. 



lower part of Fig. 2. The effective EOS appears softer 
and even more attractive compared to ground state nu- 
clear matter. This is a general feature of colhding nuclear 
matter and can be understood by a very transparent and 
model independent argument: 

Let us consider two currents with sufficiently high u, 
i.e. well separated Fermi ellipsoids in momentum space 
(01 • 82 = 0). For the following discussion we assume 
that the self-energies have no explicit momentum depen- 
dence, like in relativistic mean field (RMF) or density 
dependent RMF theory j2l] . Then the energy density (in 
the cm. frame) is given by 

T?2i9l2,u) =< E* >i2 -[^s^,gs^, +^i2gi2]/2 (10) 

and the non-locality of the system, i.e. the high relative 
momenta of the separated ellipsoids enter only via the 
momentum dependence of E*. The separation of pro- 
jectile and target nucleons in momentum space increases 
the phase space volume since in both currents states are 
occupied up to kp = For a purely local in- 

teraction which is insensitive to the relative momenta of 
the currents this increase of phase space would have the 
same effect as the introduction of an additional virtual 
degree of freedom as illustrated in Fig. 3. Thus, in a mean 
field approach the effective binding energy (| 10) can be 
approximated by a modification of the corresponding ex- 
pression for single nuclear matter 
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However, a meaningful discussion of non-equilibrium 
effects with respect to the ground state EOS should be 
based on the binding energy and thus the contribution 
from the relative motion of the two currents has to be 
subtracted. This leads to an effective EOS in colliding 
nuclear matter which is directly linked to the hydrody- 
namical picture |]l[. To do so, we subtract the kinetic 
energy of the relative motion with respect to single nu- 
clear matter at rest 



£rel{gi2,u) =< E* - M* >12 - < E* - M* >„=o 



(8) 



The kinetic energy of a nucleon inside the medium -Ejjjjj ~ 
E* - M* = Vk*2 + M*2 - M* calculated in CNM or 
NM, respectively, is thereby averaged of the correspond- 
ing configurations. Thus one obtains the binding energy 
per particle 
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as a function of the total cm. density Qyi and the cm. 
streaming velocity u. The binding energy, i.e. the ef- 
fective EOS in colliding nuclear matter, is shown in the 






FIG. 3. Schematic representation of the phase space in nu- 
clear matter (a), colliding nuclear matter (b) and in colliding 
nuclear matter as experienced by a local potential (c). 

To do so, on the right hand side of Eq. (|l^) all 
quantities, i.e. pi, S, . . . are obtained by the integration 
over one Fermi sphere at rest with the Fermi momen- 
tum kp corresponding to the density of only one current, 
however, taking thereby an enhanced phase space factor 
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K = 8 in Eq. (g). For the vector density this leads to a 
hnear dependence Qi{kp^K, = 8) = 2gi{kp,K = 4) which 
restores the total density. The dependence of the scalar 
density gsi is non-linear. However, the total cm. vec- 
tor density in colliding nuclear matter is still enhanced 
by a 7-factor, i.e. gi2{kp,kp,u) — j{u)gi{kF, k — 8) 
which does not completely cancel in the effective bind- 
ing energy and leads to a stronger repulsion originating 
form the vector field as compared to the approximation 
(pi]). In Fig. 4 the corresponding EOS ( p] ) is shown as 
obtained in the density dependent RMF approach to 
the DBHF model. Here we neglect the explicit, but weak 
p8[ momentum dependence of the DBHF self-energies. 
Varying in these calculations k from 4 to 8 illustrates the 
phase space effects. There is, of course, no exact agree- 
ment of the simplified ansatz of Eq. (O) with full self- 
consistent CNM calculations (Fig. 2), in particular at high 
densities, which is expected. With increasing density self- 
consistency effects between the effective mass M* and 
the configuration (^) itself become more important. A 
reduced in-medium mass M* leads to reduced momenta 
fc* = u^M* at fixed velocities and to a shift of the Fermi 
ellipsoids. Thus at high densities Pauli blocking effects 
persist for interacting two-Fermi-ellipsoid configurations 
even at high relative velocities . Furthermore, the rel- 
ativistic mean field approximation is non-local through 
the momentum dependence via E* in (pT[). However, it 
becomes clear from this comparison that the leading or- 
der effect for the EOS in colliding nuclear matter is the 
separation of phase space. 
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FIG. 4. EOS in the relativistic density dependent mean 
field approximation to the DBHF model. The number of de- 
grees of freedom is varied from n — 2 (neutron matter), k — 4 
(nuclear matter) to k = 6, 8. The shadowed band is limited by 
the EOSs in colliding nuclear matter with streaming velocities 
u = 0.6/0.8 (upper/lower curves). 

It is important to realize that this type of phase space 
effects is not included in standard transport calculations 
for heavy ion collisions, even when momentum dependent 
interactions are used. Phenomenological potentials M 



are usually composed by a local, density dependent po- 
tential Uioc{g) and a non-local momentum dependent 
part C/„onioc(e, k) = / d^k'f{k')V{k-k') with V an effec- 
tive two-body-interaction. In Hartree-Fock approxima- 
tion these two terms correspond to the direct (Hartree) 
and the exchange (Fock) part of the potential. In ground 
state nuclear matter the Fock terms give usuall y s mall 
corrections to the EOS (mean field dominance) [g2[. In 
colliding nuclear matter the non-local part of the interac- 
tion is responsible for the strong repulsion seen in Fig. 2. 
Applying potentials of the form ( p^ ) in transport calcu- 
lations for heavy ion collisions the Fock part J7nonioc(g, k) 
accounts by definition properly for the actual momentum 
space configurations /(k'), given e.g. by testparticle dis- 
tributions. The mean field or local part is, however, de- 
coupled from the anisotropy of the phase space since it 
is parameterized as a function of the total density. Con- 
sequently, Uiocig) reflects a density dependence which is 
only correct in equilibrated nuclear matter but does not 
apply to anisotropic momentum space configurations. 

To summarise, the equation of state probed by the 
compression phase in energetic heavy ion reactions is to 
large extent governed by local non-equilibrium. The cor- 
responding separation of phase space can be regarded as 
the introduction of an effectively new degree of freedom 
which lowers the binding energy per particle and makes 
the effective EOS seen in heavy ion reactions significantly 
softer. We conclude that this "trivial" but leading order 
effect should be taken into account when conclusions on 
the EOS arc drawn from heavy ion collisions. 
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